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Abstract 



A self-contained Fortran-90 program based on a classical trajectory model with 
stochastic breakup is presented, which should be a powerful tool for quantifying 
complete and incomplete fusion, and breakup in reactions induced by weakly-bound 
two-body projectiles near the Coulomb barrier. The code calculates complete and 
incomplete fusion cross sections and their angular momentum distribution, as well 
^ ■ as breakup observables (angle, kinetic energy and relative energy distributions). 
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cross section, Spin distribution, Kinetic energy distribution, Relative energy distri- 
bution, Angle distribution, Classical trajectory 
PACS: 25.70.Jj, 25.70.Mn 

Classification: 

External routines /libraries: 

Several source routines from Numerical Recipies, and the Mersenne Twister random 
number generator package are included to enable independent compilation. 
Nature of problem: 

The program calculates complete and incomplete fusion cross sections and their spin 
distribution, as well as breakup observables (e.g. the angle, kinetic energy, and rel- 
ative energy distributions of the fragments) in reactions induced by weakly-bound 
two-body nuclei near the Coulomb barrier. 
Solution method: 

A classical trajectory model with stochastic breakup is used to calculate all the 

observables. See Ref. [1] for further details. 

Restrictions: 

The program is suited for weakly-bound two-body projectiles. The initial orienta- 
tion of the segment joining the two breakup fragments is considered to be isotropic. 
Running time: 

About one hour for input provided, using a PC with 1.5 GHz processor. 
References: 

[1] A. Diaz-Torres et al, Phys. Rev. Lett. 98(2007) 152701. 



1 Introduction 

Recent developments of radioactive isotope accelerators enable the investiga- 
tion of fusion reactions that form heavy elements in the cosmos. These involve 
reactions of nuclei far from stability. The most exotic of these are often very 
weakly bound. Breakup of weakly bound nuclei is an important process in 
their reactions with other nuclei. A major consequence of breakup is that not 
all the resulting breakup fragments might be captured by the target, termed 
incomplete fusion (icf); capture of the entire projectile by the target is called 
complete fusion (cf). Such icf processes can dramatically change the na- 
ture of the reaction products, as has been investigated in detail for the stable 
weakly-bound nuclei 9 Be and 6 ' 7 Li [I]. There, at energies above the fusion bar- 
rier, CF yields were found to be only ~ 2/3 of those expected, the remaining 
1/3 being in icf products. Events where the projectile breaks up and none of 
the fragments are captured provide an important diagnostic of the reaction 
dynamics. This we call no-capture breakup (ncbu), which is expected to be 
predominant at energies below the fusion barrier. The modelling of all these 
reaction processes within the same theoretical framework is an outstanding 
problem. The continuum-discretised-coupled-channels (CDCC) quantum me- 
chanical model can make reliable predictions of the ncbu process, but cannot 
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distinguish between icf and CF processes unambiguously |2J. An alternative 
for solving this problem is the development of classical dynamical approaches 
based on the concept of classical trajectories, which are being successfully 
applied to describe atomic many-body collisions, including ionization [1], a 
process in some ways analogous to breakup. 

In this paper, we present a program of the classical trajectory model with 
stochastic breakup, recently published by the author et al. in Ref. [3]. Since 
the code could be very useful for many researchers involved in near-barrier 
experiments of fusion and breakup with weakly-bound/ radioactive nuclei, it 
will be beneficial make the program accessible to everyone. In Section 2, a 
description of the model is given to make the paper self-contained, whilst in 
Section 3 the program and the input file are explained. The code is illustrated 
in section 4 with the reaction of a pseudo- 8 Be projectile P (assuming a weakly- 
bound state of two a-particles) with a 208 Pb target. For this test case, the good 
agreement of the classical model calculations for the ncbu process with those 
of the CDCC quantum mechanical model shows the reliability of the classical 
approach (see Ref. [3] for details). 



2 Model 

2.1 Projectile-target interaction 

The weakly bound (two-body) projectile P, with incident energy E and or- 
bital angular momentum L , is incident on the target T, initially at rest in 
the laboratory frame. Prior to breakup the projectile follows a classical orbit, 
which is characterised by a P-T distance of closest approach R min (E , L ). 
These orbits of the bound projectile are calculated by numerical solution of 
the Newtonian equations of motion in the presence of the Coulomb and nu- 
clear potentials between P and T. These generate a fusion barrier for head-on 
(Lq = 0) collisions having a height V^ T at an internuclear distance R^ T - 

2. 2 Encoding of breakup 

The Coulomb and nuclear two-body interactions that cause breakup are em- 
pirically encoded in a density of (local breakup) probability Vbu(R), which is 
a function only of the projectile-target separation R. It is defined such that 
V^u{R)dR is the probability of breakup in the interval R to R + dR. A key 
feature that must be emphasised is that for a given projectile-target combi- 
nation, both experimental data [5] and CDCC calculations indicate that the 
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integral of this breakup probability along a given classical orbit is an expo- 
nential function of its distance of closest approach, R min (E , L ). Explicitly 
(see Appendix A), 



oc 




(1) 



-ft? 



It follows uniquely from this equation that the local breakup function has the 
same exponential form, 



2. 3 Initial conditions of breakup events 

The code then proceeds with a Monte Carlo approach to sample the initial con- 
ditions of breakup. For each L (chosen to be an integer number of K) a sample 
of N incident projectiles is taken. The position of breakup on this orbit is de- 
termined by sampling a breakup radius Rbu in the interval [R m i n (E , L ), oo] 
with the weighting V^ U (R). Although dependent on the constants A and a, 
the exponential weighting of Vbu(R) will clearly place most Rbu m the vicin- 
ity of R m in- If the chosen L is less than the critical partial wave for projectile 
fusion, L cr , then the associated trajectory would normally lead to CF, i.e. 
Rmin < Rb T - F° r these Lq we set R m in = Rb T , when sampling Rbu, and all 
breakup events are confined to the incoming branch of the projectile trajec- 
tory. For L > L cr breakup can take place on both the entrance and the exit 
branches of the classical orbit, which are sampled equally. Having chosen the 
position with the projectile center of mass at breakup, it is instantaneously 
broken up into fragments Fl and F2. 

Following breakup, the two fragments Fl and F2 interact with T, and with 
each other, through real central two-body potentials having Coulomb barri- 
ers V B 3 at separations Rg, i,j = 1,2,T, % ^ j. The instantaneous dynamical 
variables of the excited projectile at breakup, namely its total internal energy 
£12, its angular momentum £ 12 and the separation of the fragments d\ 2 are all 
Monte Carlo sampled. The initial separations di 2 between the fragments are 
Gaussian distributed in their classically allowed region. This mimics the radial 
probability distribution of the projectile ground-state (g.s.) wave function. For 
high £ 12 , when there is no barrier between Fl and F2, d 12 is taken as their 
external turning point. In the test calculation described below, the orientation 

— * — * 

of d\2 is chosen randomly over the 4ir solid angle, and the orientation of £12 
is chosen randomly from all directions orthogonal to d\ 2 - Other distributions 
could be taken in cases where the two breakup fragments are not identical. t\ 2 
is sampled uniformly in the interval [0,£ ma:r ]. For e 12 two sampling functions 




(2) 



4 



were tested ranging from the energy of the top of the barrier, Vg 2 , to a cho- 
sen maximum e max . Using uniform weighting or an exponentially decreasing 
weighting yielded very similar outcomes. The convergence of the observables 
was faster for the latter case, which was therefore used. Both £ max and e max 
were increased until convergence of the observables occured, as in the CDCC 
calculations. 

2.4 Trajectories of breakup fragments and the target 

Having fixed the positions and dynamical variables of the excited projectile 
fragments at the moment of breakup, the instantaneous velocities of the parti- 
cles Fl, F2 and T are determined by conservation of energy, linear momentum 
and angular momentum in the overall center of mass frame of the projectile 
and target system (see Appendix B for details). These breakup initial condi- 
tions were transformed to the laboratory frame where the equations of motion 
are solved. The calculated trajectories of Fl, F2 and T determine the number 
of ICF, CF and ncbu events, fragment Fj being assumed to be captured if the 
classical trajectories take it within the fragment-target barrier radius R^ . 

2.5 Probabilities and cross sections 

From the N breakup events sampled for each projectile angular momentum L , 
the numbers of events N { in which i = (ncbu), 1 (icf), or 2 (cf) fragments 
are captured determine the relative yields Pi = Ni/N of these three reaction 
processes after breakup, with P + Pi + P 2 = 1. The absolute probabilities 
Pi(E , Lq) of these processes are expressed in terms of the relative yields and 
the integrated breakup probability over the whole trajectory PBu(Rmin)'- 



where H(x) is the Heaviside step function. The cross sections are calculated 
using 



Po(E ,L ) 

Pi{E ,L ) 
P2{E ,L ) 



— PBu(Rmin) Pq, 
= PBu(Rmin) Pi, 




(3) 
(4) 



(5) 



<Ti(E ) = vrA 2 ^(2L + 1)^(^0,^0), 



(6) 



where A 2 = H 2 



/[2mpE ] and mp is the projectile mass. 
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The other observables, such as the angle, kinetic energy and relative energy 
distributions of the fragments from ncbu events, are calculated by tracking 
their trajectories to a large distance from the target, which was 200 fm in the 
calculations given below. 



3 Computer program and input file 

3. 1 Structure of the code 

The code has a main program and eight modules. The main program breakupSD 
directs the input to be read, the problem to be solved, and details of the cal- 
culation to be written in output files. The modules are global data, potentials, 
mtl9937, nrutil, initial conditions, fusion, Angular momentum distribution 
and INPUT VALUES. 

The main program breakupSD calls the module INPUT VALUES first, in 
which the subroutine input data reads the input file described below. For a 
given partial wave between the projectile and the target (IMPACTMIN up to 
IMPACTMAX), the subroutines projectile trajectory and trajectory arrays of 
this module calculate the orbit of the bound projectile and store it for inter- 
polations in the module initial conditions. Thereafter, ISEEDMAX breakup 
events with sampled initial conditions are calculated for every partial wave. 
The initial conditions for the propagation in time of the three bodies are fixed 
by the subroutine initial values of the module initial conditions, whilst the 
classical trajectory of the breakup fragments and the target are calculated by 
the driver-subroutine ODEINT of the module nrutil. During the propagation, 
the possible capture of the breakup fragments by the target is determined 
by the subroutine fusion events in the module fusion. Finally, the spin dis- 
tribution and cross sections for CF and ICF, and for the breakup process are 
calculated by the subroutine spin distribution in the module Angular momen- 
tum distribution. Here, other observables such as the angle, kinetic energy and 
relative energy distributions of the fragments from ncbu events are also com- 
puted. Depending on the value of the output control variables (FILE1, FILE2, 
FILE3 and FILE4), details of the calculation can be written into output files. 

Module global data. This defines global variables used in different modules. 

Module potentials. Here the nuclear and Coulomb interactions between the 
participants of the reaction (projectile fragments and the target) are defined. 

Module mtl9937. It contains the Mersenne Twister random number generator, 
written by Makoto Matsumoto and Takuji Nishimura [BJ. 
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Module nrutil. In this module several subroutines from Numerical Recipies 
are included, which are mainly related to integrating the classical equations 
of motion. ODEINT is the integrator driver that calls the subroutine RKQC 
which is a fouth-order Runge-Kutta integrator ensuring accuracy and adjusted 
stepsize. The forces are defined in the subroutine DERIVS. The module also 
contains subroutines (printout 1 and printout2) for writing into output files 
the trajectory of the nuclei along with the value of the integrals of motion. 
The latter is crucial for checking the accuracy of the integration of the classical 
equations of motion. 

Module initial conditions. This module defines the initial conditions for the 
propagation in time of the three bodies, following the breakup of the two- 
body projectile. This task is performed by the subroutine initial values, which 
is guided by the conservation of energy, linear momentum and angular momen- 
tum in the overall center-of-mass frame of the projectile and target system. 
This is transformed to the laboratory reference frame (with Galilean kine- 
matical relations) where the equations of motion are solved using a system 
of spherical coordinates. This module calls the module mtl9937 in the sam- 
pling of (i) the initial excitation energy of the projectile, (ii) the initial relative 
angular momentum between the projectile fragments, (iii) the initial separa- 
tion between the fragments, (iv) the breakup radius, (v) the initial orientation 
of the segment joining the two fragments of the projectile, and (vi) the ini- 
tial direction of the radial velocity along that segment. If the output control 
variables are activated, the initial conditions will be written into output files. 

Module fusion. In this module the subroutine fusion events analyses the pos- 
sible capture of any of the breakup fragments by the target. Here, the relative 
energy and angular momentum between the nuclei during the propagation are 
calculated by the subroutines relative energy and relative spin. 

Module Angular momentum distribution. Knowing the statistics of fusion and 
ncbu events after a large number of sampling (1000 per partial wave in the 
example given below), the subroutine spin distribution calculates observables 
related to CF and ICF processes (spin distribution and cross sections). The 
angle, kinetic energy and relative energy distributions of the ncbu events and 
the breakup cross section are also computed here. This subroutine extensively 
calls the subroutine LOCATE of the module nrutil. All observables are written 
into output files. 

Module INPUT VALUES. Here the input file is read by the subroutine input 
data. It also calculates the s-barrier features (radius and height) between (i) 
the projectile and the target, (ii) the projectile fragments and the target, 
and (iii) the two fragments of the projectile. The module also contains the 
subroutines projectile trajectory and trajectory arrays mentioned above. 
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3.2 Input file 



Fig. [T] shows the input file of the program. At the bottom of the figure, the 
namelist of the input variables appearing in the code are shown. The lines 
related to the potentials (lines 12-19) are self-explanatory, whilst some vari- 
ables in previous lines have already been mentioned in the description of the 
program. We will here describe lines 1-10 only. 

The integer variables of the first line allow the user to write details of the 
dynamical calculations into output files, when they are equal to one. FILE1 
and FILE2 open output files to write the trajectory of the breakup fragments 
and the initial projectile-target, respectively. The file opened by FILE3 con- 
tains details of breakup events, and FILE4 is for plotting the trajectory of the 
breakup fragments. Line 2 defines the window of orbital angular momentum of 
the incident projectile (IMPACTMIN, IMPACTMAX) in units of h, and the 
number of breakup events per partial wave (ISEEDMAX). PARTICULAR 
IMPACT in line 3 is to select a partial wave, whose associated dynamical 
calculations are written into the file opened by FILE3. Line 4 refers to the 
incident energy of the projectile in MeV (E0). Lines 5-6 define the range of 
initial excitation energy (EXCMIN, EXCMAX) in MeV and relative angu- 
lar momentum (L12MIN, L12MAX) in units of h for the breakup fragments. 
In line 7, TYPE EXC controls the sampling function for excitation where for 
TYPE EXC=0 the weighting is uniform, and for TYPE EXC=1 it is exponen- 
tially decreasing with a coefficient alfaexc. Line 8 defines the centroid (d012) 
and width (sig012) in fm of the Gaussian that describes the radial probability 
distribution of projectile ground-state wave function. It is used to sample the 
initial separation between the breakup fragments. The breakup function [see 
expression (TjQ) in Section 2] is given in line 9 by the parameters alpha and beta, 
being beta=ln(A). Line 10 defines a maximal projectile-target separation for 
sampling the breakup radius. 



4 Test run 

Using the input file of Fig. [H Figs. El [3] and 0] show output files for angular mo- 
mentum distribution in icf, CF and ncbu processes, respectively. The angular 
momenta are in the first column, whilst in the second and third columns are 
partial probabilities and cross sections. The total cross section appears at the 
end of the file. In Fig. [21 the ICF spin distribution refers to either the angular 
momentum brought by the captured fragment into the target or the orbital 
angular momentum of the incident projectile. The former has been employed 
very recently to understand isomer ratio measurements [7J. Of course, both 
distributions are the same for CF. Since the component of the CF spin distri- 
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<m <U5 otono 

000 
(J50.D0 
(fc.UO urn 

(Sl.fi5 DO.O0 

1 m.922 
oi. nii OigC 

0.9220 9.73O0 
050.00 

«8fttt8##POTENTIAlJi## tf^JlillJ)llJIJIJlUJIIIHttt 
20H.U 082.0 
(XJH.il QO4.0 
(XH.O 0O2.0 
(XH.fl 002.0 

-120.WGO li'SjQO (S.7552 U2O00 

-062.UTO I.390CJ (1.6200 1:2000 

JESSZ.CWW 1.3900 O.fi200 I. HNS) 

#Jjt,S»64 1 .2000 fltfiSO 1.2000 

LiBS J . F1LEJ, F1LE2, F1LE3, F1LE4 

Line 2: [MPACTMIN, JMPACTMAX. ISEEDMAX 

Line 3: PARTICULAR IMPACT 

Line 4: Efl 

Line 5: EXCMAX, L12MAX 
Line ft EXCMlhi. Lt2MIN 

Line 7: TYPE, EXC, aKasjie l rally impi-mam wheu TYPE,EXC=J f 
Line S;dOI2.aigHI2 
Line 9: alpha. beLs 
Line J ft RBU , max 

.Line 12: MaMi and eliacee of (he emj^i 1 AT.ZT) 
Line LI: Ma^ ;uiJ lIlu'uu oriliE jwsjjceiLL' ( AP.Zf i 
Line U: Mass, and: ehargc of ihe. fngmejvl I fAP J JZPJ } 
Line 15: Mas.ii in J Llurce <s( lhe fragment (A.P2,ZP2) 
Line 16. 'J'itf^i-J'iujLviilt- fKHcriijiil iVdTP.jrriJTP..iDTPi 

and Cnukimi) radius f ireflTP) 
Line 1 7. Turgtci-Fi ugitienL L fhncjiuul iVOI.irill.jdli and 

and Coutomb radium frreOt i 
Line JH: Tflrgct-Fragmtnil (KMejuiiU lV(E,frfl2.afJ2i and 

and Coulomb radius i;rre02 1 
Line J a* Fragmcntf- Fragment potential ( VOL2.nO]2,aE3l2l iad 

and Coulomb radius rrre0t2i 

Fig. 1. Input file for PLATYPUS code (PL ATYPUSinp) . See text for further details. 

bution associated with the first term of expression (jSJ) (fusion of the bound 
projectile) is clearly equal to 7rA 2 (2L + 1)[1 — PBu(Rmin)], in Fig.[3]we present 
only the nontrivial component related to the capture of all the projectile frag- 
ments after breakup. The code also calculates CF cross sections related to the 
inert projectile and to the mentioned two components of this process when 
the projectile can be dissociated, as shown at the end of Fig. [3j As expected, 
the ncbu process in Fig. H] reveals a broader angular momentum distribution 
than the CF and ICF processes. 

A Appendix 
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REGARDING ANGULAR MOMENTUM BROUGHT 
BY FRAGMENTS rNTO TARGET 



LXlbar) P_U SIGMA2_L (mb) 

****** SMOOTH DISTRIBUTIONS***** 



0.00 .725E-02 0.S7O5EJ01 

I. 00 .S57E-01 0.1232E-HH 
2.00 .221E-KM 0.4141E4fli 
3.00 .423E+00 D.9467E-H31 
4.00 .635E+00 0.1611E+02 
5.00 .756E+O0 0.2147E4O2 
6.00 .M2E+00 0.3059E-H32 
7.00 .102E-KI1 0.3729E+02 
8.00 .9S3E+00 0.41JSE-H32 
9.00 762E+00 D.3425E-H32 
10.00 .454E-0C 0.2106E-KJ2 

II. 00 .274E-00 Q.1264E-+32 
12.00 .177E-00 0.7912E+01 
13.00 .100E-00 0.453 3E-KU 
14.00 .357E-Q1 0.1594E+01 
15.00 .178E-01 0.SO71E-+OQ 
16.00 .396E-02 0.187OE-+AQ 
17.00 .132E-02 0.6663E-O1 



TOTAL SIGMA ICT =0.2446E-03 mb 



REGARDrNG PROJECTILE ANGULAR MOMENTUM 

0.S439E-O1 0.1376E4O0 

1 0.:<?96E-K)0 0.1319E+01 

2 0.:657E-H)0 0.21&6E+01 

3 0.3019E-KJ0 0.3446E-KJ1 

4 0^70E-H)0 0.3M7E4O1 

5 0.3039E-K10 0.:450E-K)1 

6 0.3072E-KJ0 0.6511E+01 

7 0.33S5E-K10 0.83O2E-K)l 

8 0.3474E-HJ0 Q.9G29E401 

9 0.4134E-KJ0 0.12BOE-H52 

10 0.4127E+00 0.1413E+O2 

11 0.4200E-0G 0.1573E4O2 

12 0.4740E-00 Q.1932E-KI2 

13 0.4E13E-00 0.211SE4O2 

14 0.5070E-OO 0.2397E-KJ2 

15 0.45E2E-00 0.2316E+O2 

16 0.3E96E-00 0.2096E-+O2 

17 0.3323E-00 0.1896E-KI2 

18 0.2657E-0G 0.16D3E-KT2 

19 0.1160E-00 Q.7377E-KU 

20 0.570SE-01 0.3S15E+01 

21 0.4295E-01 0.30UE+01 

22 0.197QE-Q1 Q.1445E-KH 

23 0.137EE-01 0.1056E+O1 

24 0.5Q48E-02 0.403 3E-+O0 

25 0.25GQE-Q2 0.2D7BE-+O0 

26 0.6901E-03 0.5963E-O1 

27 0.45S2E-03 0.410BE-01 

28 0.1903E-03 0.1769E-O1 

29 0.1582E-03 0.1521EM>1 

SIGMA2_ICF=0.2446E+03 mb 



Fig. 2. icf angular momentum distribution (ICF SPIN DISTRIBUTION). See text 
for further details. 



REGARDING AKGLIAR MOMENTUM BROUGHT 
BY FRAGMENTS INTO TARGET 

L(bbar) P_L SIGMA2J. (mb) 



****** SMOOTH DISTRIBUTIONS***** 



0.00 .572E+O0 0.931BE-+O0 
LOO .3B5E+00 0.1SS3E-KI1 
2.00 .3S8E+O0 0.316OE+O1 
3.00 .354E+00 D.4040E-+A1 
4.00 .384E+00 0.563OE-KU 
5.00 .353E+O0 D.6325E401 
6.00 344R+0Q 0.7293E-KH 
7.00 J12E4O0 0.7641E-KH 
8.00 .299E+QQ 0.8295E4O1 
9.00 .230E+QO 0.7127E-+O1 
10. OS .223E-00 0.765 1E-KH 
11.00 .20SE-00 0.781 1E+01 
12.00 .142E-00 0.5804E-H01 
13.00 .101E-00 0.4440E-KI1 
14.00 .475E-01 0.2244E-+O1 
15.00 .21SE-01 0.1O99E-+O1 
16.00 .165E-01 D.SS67E4O0 
17.00 .557E-02 D.489CE-H30 
IS. 00 .593E-02 0.3579E4A0 
19.00 .19SE-02 0.1257E-W0 



*** *** **" t* 1 1* **» M * M «tt *** t* i** t 



TOTAL SIGMA CF =0.8323E-D2 mb 



5IGMA_CF (INERT PRO J>=0 .652 1E-0 3 mb 
SIGMA2_CF=0.£323E-02 mb 
SIGMA_CFCNO BU)=0.2222E+03 mb 
SIGMA_CF(TOTH).3054E+03 mb 



Fig. 3. CF angular momentum distribution (CF SPIN DISTRIBUTION). See text 
for further details. 

Let us define two probabilities: (i) the probability of breakup between R and 
R + dR, p(R)dR [being p(R) a density of probability], and (ii) the probability 
the weakly-bound projectile has survived from oo to R, S(R). The survival 
probability at R + dR, S(R + dR), can be written as follows 

S(R + dR) = S(R) [1 - p(R)dR] . (7) 



Expression (JTj) suggests the following differential equation for the survival 
probability S(R), 

= -S{R) p(R), (8) 
whose solution is [S'(oo) = 1]: 

R 

S(R) = exp(- J p(R)dR). (9) 



11 



REGARDING PROJECTILE ANGULAR MOMENTUM- 



L(hbar) PBU_L SIGMAL 

0.3:9 GE-02 0.5374E-D2 

1 0.461 5&J02 G.2257E-01 

2 0.5934E-02 Q.4S36E-D1 

3 0.3296E-O2 0.3 762E-01 

4 8571E-Q2 O.1257E-KJ0 
: 0.2637E-G2 0.472 9E-D1 

6 0.791 l&W 0.1677E-K30 

7 0.7252E-Q2 0.1773E-K3Q 
E 0.1253E-O1 0.3471E-+O0 

9 0.1582E-O1 0.490 1E+0Q 

10 0.2307E-01 0.7S99E+D0 

11 D.3099E-01 0.1162E4A1 

12 D.4255E-01 0.1746E-H31 

13 D.7714E-01 0.33?5E+01 

14 D.1048E-0C 0.4956E-H31 

15 0.1793E-00 G.9062E-KH 

16 D.2532E-DS 0.1362E-K)2 

17 0.31S4E-0G 0.1817E+O2 
IS 0.3S77E-00 Q.233SE-K32 

19 0.i327E-00 0.3387E-+O2 

20 D.4254E-D0 0.2843E-H32 

21 D.331SE-DO Q.2326E-K32 

22 0.2740E-00 Q.2010E-HJ2 

23 0.22S7E-00 0.1752E-KJ2 

24 0.1930E+00 Q.1542E-K32 

25 D.159QE-D0 0.1322E+02 

26 01337E+00 Q.1156E-KJ2 

27 0.1115E-00 Q.9993E-KH 

28 D.9050E-Q1 O.S410E+01 

29 D.7569E-01 0.728OE-+A1 

30 D.6349E-01 0.6314E-KH 

31 0.5241E-01 0.53S3E-+O1 

32 0.4346E-01 0.4606E-+O1 

33 D.362SE-01 0.3963E+D1 

34 D.2994E-01 0.336SE+D1 

35 0.2495E-01 0.2885E-+A1 

36 D.2062E-01 0.2453E-KU 

37 0.1704E-01 0.20S3E-K31 

38 0.1407E^1 0.1766E-KI1 

39 0.1157E-01 0.1490E-KH 

40 0.955 1E-42 0.12S1E+01 

41 D.7££0E-O2 0.1D6&E+31 

42 D.M82E^2 O.S9S2E-K)Q 

43 0.5350E-02 0.75S7E+O0 

44 D.439SE-02 0.6381E-KIQ 

45 D.3553E-02 0.5271E-+flO 



SIGMA 3U=:.3063E-03ii± 



Fig. 4. ncbu angular momentum distribution (BU SPIN DISTRIBUTION), 
text for further details. 
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From ©, the breakup probability at R,B(R) = 1- S(R). If J*p(R)dR <C 1, 
B{R) can be written as 

R 

B(R) « y p(R)dR. (10) 



From (flQj) . identifying p(-R) with V^ U (R), we obtain expression (TjQ) for the 
breakup probability integrated along a given classical orbit. 

B Appendix 

The integrals of motion in the overall center-of-mass (CM) system are the 
total energy E tot = j^^^Eq, the total linear momentum P tot = 0, and the 

total angular momentum L tot = mpb (v — Vcm) that is orthogonal to the 
initial reaction plane, tut, vnp, bo, v, and Vcm are the mass of the target 
and projectile, the impact parameter between the projectile and the target, 
the velocity of the incident projectile in the laboratory system and the CM 
velocity, respectively. 

Just after breakup, the two-body projectile is excited to a definite state (£% 2 , 
£12 and d.12), as explained in Section 2. The relative vector between P and T 
(Rpt) is also known. Thus, the separation between the three bodies is known. 
The modulus of the velocity between P and T (Vpt = Ppt/ppt) results from 
the total energy conservation 

E tot = e 12 + Urr(r 1T ) + U 2T {r 2 T) + Pp T /2/j PTj (11) 



where U is the interaction potential between the target and the breakup frag- 
ments. 

The total linear momentum P tot = Pt + Pi + P2 — Pt + Pp*, where pp* is 
the momentum of the center of mass of excited P relative to the overall CM. 
We need the velocities of P and T relative to the overall CM (vp and Vj) to 
complete the initial conditions for subsequent propagation in time of the three 
bodies. These velocities are related to each other by the expressions 

mp ^ 

v T = v P , (12) 

rriT 

Vpt = Vp — Vt, (13) 



where the magnitude of Vpt is known through expression (fTTl) . To know the 
direction of this velocity the conservation of total angular momentum is ap- 
plied. 
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The total angular momentum L tot = £i 2 + Lpp, so the angular momentum 
(Lpt) associated with relative motion of P and T about CM is known. This 
vector can be written as 

Lp T = mpRpT x vp. (14) 

We now write vp in terms of radial and transverse components as follows: 
Vp = v^p r + v^p q, (15) 



where r = R PT /R PT and q = n x r, being n = L PT /L PT . The transverse 
component v^p = L PT /(mpRpT), and for the target w§' = —Lpt/(ititRpt)- 
The radial component is obtained using expressions ( fT2l) -( fT3l) and knowing 
the transverse component: 



~(r) 

v P 



±{V* T - 



Vp 



i+— )l 2 } 1/2 /(i+— ) 



(16) 



Both positive and negative roots are consistent with the conservation of the 
integrals of motion. Hence, both roots are uniformly sampled. Finally, the 
position and velocity vectors of the projectile fragments and the target are 
transformed to the laboratory system using Galilean transformations. 
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